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ABSTRACT 


An analysis is presented of low-level trajectory data obtained by means of constant level balloon flights from 
Cape Hatteras, N.C., during September and October 1959. An approximately constant floating level was obtained 
by flying the nearly constant volume Mylar balloons (tetroons) with an internal superpressure of about 100 mb. 
The metalized tetroons were positioned at l-min. intervals by means of a manually operated SP-1M radar. From 
knowledge of these positions, overlapping 5-min. average velocities were determined for flight durations of up to 
5 hr. On some of the flights the radar return was enhanced by the addition of a radar reflective mesh to the tetroon. 
With the addition of this mesh, flights at altitudes of less than 5,000 ft. were tracked as far as 92 n. mi. from the radar, 
or approximately to the radar horizon. 

Spectral analysis of the velocity data obtained from the four best flights shows some evidence for a (Lagrangian) 
wind speed periodicity of 26-min. period, a vertical motion periodicity of 13-min. period, and a cross-stream velocity 
periodicity of 17-min. period. Cross spectrum analysis shows that, with the exception of oscillations of 45-min. period, 
the wind speed is at a maximum ahead of the trough in the trajectory. Thus, if the large-scale air flow is nearly 
geostrophic, there is evidence that kinetic energy and momentum are transported down the pressure gradient by 
these small-scale oscillations. The maximum upward motion of the tetroon tends to take place near the trajectory 
trough line for oscillations of a period exceeding 18 min. and near the trajectory crest for oscillations of smaller period. 
Therefore, looking downstream, the longer-period tetroon oscillations tend to be counterclockwise in a plane normal 
to the mean trajectory while the shorter-period oscillations tend to be clockwise. However, until more information 
is obtained on the small-scale temperature field and the extent to which the tetroons follow the vertical air motion, 
any statement regarding ‘“‘direct’’ and “‘indirect”’ air circulations is tentative. 

The ratio of one minus the cross-stream and one minus the along-stream autocorrelation coefficients for these 
flights is approximately 0.6, suggesting a certain similarity between Eulerian space and Lagrangian autocorrelation 
coefficients. The tetroon data also indicate that initially one minus the cross-stream autocorrelation coefficient is 
proportional to the time, as would be anticipated from Lagrangian turbulence theory. These data tend to confirm 
that, for the scale of motion under consideration, the Lagrangian-Eulerian scale factor 8 of Hay and Pasquill has 
a value near 4, 
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1. INTRODUCTION the usefulness of such balloons for the delineation of the 


large-scale air flow at jet stream levels over a large seg- 


i of the United States Navy has demonstrated 


formed in connection with Weather Bureau research for the U.S. Atomic 


ment of the Northern Hemisphere [{1, 2]. The purpose of 
this article is to show how constant level balloons may be 
utilized for the delineation of the small-scale air flow at 
levels near the surface and at distances up to 100 n. mi. 


mission and the U.S, Public Health Servies. from a radar site. The Lagrangian fluctuations derived 
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from such low-level flights are of great interest in turbu- 
lence and diffusion studies. 


2. THE SUPERPRESSURED BALLOON CONCEPT 


Previous attempts to obtain the Lagrangian charac- 
teristics of the air flow at low levels have, in most cases, 
been based on the use of “neutral” or “no-lift’’ elastic 
balloons. While information of interest has been ob- 
tained from such flights [7], they have not proved particu- 
larly suitable for the determination of the Lagrangian 
characteristics of the flow over relatively long distances. 
Since constant volume, superpressured balloons have the 
capability of floating along a constant density surface 
almost indefinitely (barring large vertical air motions), it 
is appropriate to consider their use for this problem. 

A balloon ascends or descends in the atmosphere accord- 
ing to whether its bouyancy force, V,(p,—p,), (where V, 
is balloon volume, p, is air density, and p, is the density of 
a gas lighter than air, in this instance, helium) exceeds or is 
less than the weight, W>, of the balloon system. However, 
if the balloon is of constant volume then, for a given mass 
of helium introduced into the balloon at the earth’s surface, 
there is a certain air density where the bouyancy force 
and the weight of the balloon system balance. The 3- 
dimensional surface along which the atmosphere possesses 
this density is the surface along which the constant vol- 
ume balloon will float. The balloon will continue to float 
along this surface as long as it remains at full volume. The 
balloon may go slack owing to seepage of helium through 
the skin of the balloon, through nighttime cooling of the 
balloon and inclosed helium due to radiation fluxes, or 
through cooling of the balloon due to any other causes. 
For these reasons it is desirable that initially, at flight 
level, the helium within the balloon exert a greater pres- 
sure than that of the ambient air. This pressure excess 
is known as superpressure. Analytically, the helium 
pressure (p,) required to prevent a slack balloon at the 
cold temperature 7, is given by 


Pr= (T/T) Da (1) 


where p, is the ambient air pressure and 7), is the warm 
temperature. For example, if one assumes that a con- 
stant volume balloon flying at 900 mb. possesses a tem- 
perature of 300° A. during the day and a temperature 
of 280° A. during the night, then a daytime superpressure 
of 63 mb. would suffice to keep the balloon fully inflated 
at night. While statistics on day-night balloon tem- 
peratures at low flight levels are not vet available, one 
would expect that a daytime superpressure of 100 mb. 
would be suffiicient to carry the balloon through the night- 
time hours with no deviation from the density surface. 


3. EQUIPMENT AND FLIGHT TECHNIQUES 
A relatively new polyester film called Mylar ? is highly 


2 Mylar—a trade name for a polyester film manufactured by the E. I. du Pont Com- 
pany of Wilmington, Del. 
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suitable for use in the construction of constant volume 
balloons since it possesses very low permeability and high 
tensile strength (13,000 lb. in.~*). It can be shown that 
a spherical Mylar balloon of 1-foot radius and a skin 
thickness of only 2 mils (2/1000 in.) is capable of support- 
ing a superpressure of 150 mb., thus undoubtedly satis. 
fying the day-night requirements mentioned aboye. 
Furthermore, while it was originally feared that the 
seals or gores on the balloon would be unable to withstand 
such a superpressure, it was found that the G. T. Schjeldahl 
Company of Northfield, Minn. had perfected a method 
of heat sealing which yielded a seal strength comparable 
to that of the Mylar itself. For purposes of economy 
and reliability, however, it was desirable that these seals 
be straight lines. Therefore, tetrahedron-shaped balloons 
rather than spherical balloons were constructed for our 
experiments. Combination of the word “tetrahedron” 
and the word “balloon” results in the word “‘tetroon” 
as a designator for these particular balloons. 

Tetroons with a side length of 36, 42, and 60 in. (nomi- 
nal volumes, respectively, of 0.20, 0.32, 0.94 m.*) have 
been used in our experiments. It is desirable that the 
tetroon volume not exceed about 3 m.° partially in order 
to comply with regulations of the Federal Aviation 
Agency concerning the flying of balloons within air space. 
With no weight attached, the ceiling of the 42-in. tetroon 
is about 10,000 ft. while the ceiling of the 60-in. tetroon 
is about 20,000 ft. Figure 1 shows a picture of a 42-in. 
tetroon with inflation rig and manometer (for measuring 
superpressure) attached. Note that the tetroon has 
been coated with a molecular film of aluminum in order 
to provide a radar target. 

Some of the earlier tetroon flights suggested that the 
tetroon volume was not being conserved under condi- 
tions of considerable superpressure, since the tetroons 
were flying at a higher elevation than would be antici- 
pated from the balance of bouyancy force and weight on 
the assumption of constant balloon volume. If the tet- 
roon volume is not absolutely constant as it ascends 
through the atmosphere it can be seen, by equating the 
weight of the balloon system and the bouyancy force, 
that the weight, X, which must be attached to the tetroon 
to enable it to fly at a given density surface is 


X= V,(ea— pn) — Wo= V pa— Virns— Wo (2) 


where V, is the tetroon volume at the earth’s surface, 
V.=V,+AV is the tetroon volume at flight level, p, is 
the air density at flight level, p, is the helium density 
at flight level, p,, is the helium density at the earth’s 
surface, and W’, is the weight of the balloon system. 

Differentiating equation (2) with respect to pressure it 
is found that 


) ee ae A (3) 


Sp" “sp t bp 


The right hand term in this equation is negative since 
| £ 
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tetroon volume increases with decrease in ambient pres- 
sure (increase in superpressure). If this negative term 
becomes larger in magnitude than the middle term in 
equation (3), the equation states that more weight must 
be attached to the tetroon to make it float at a low pres- 
sure (high level) than at a high pressure (low level). 
Under these conditions the tetroon is not in stable equi- 
librium at any level but continuously ascends, thus in no 
sense performing the function of a constant level balloon. 
This takes place when the volume increase of the tetroon 
associated with increase in superpressure more than offsets 
the decrease in bouyancy force brought about by the 
decrease, with decrease in ambient air pressure, of the 
density difference between helium and air. It is apparent 
that the change of tetroon volume with increase in super- 
pressure must be determined precisely, not only to permit 
the flying of the tetroons at the required altitude, but 
also to ensure that a constant level balloon flight actually 
resul 

The change in tetroon volume as a function of super- 
pressure was determined by careful weigh-off at the 
ground. The amount of superpressure within the tetroon 
Was tieasured by a sensitive manometer while the tetroon 


Photograph of 42-in. tetroon with inflation rig and manometer attached. 


volume, V,, was determined from the equation 


, Wer 
V, I (4) 
Pas Phs 


where L is the free lift of the tetroon (measured by sen- 
sitive scales), p,, is surface air density, and the other 
Each of the 


tetroons tested was taken in steps up to a superpressure 


parameters are as defined previously. 
of 150 mb. with the tetroon volume being determined at 
each step. The tetroons were then deflated (a vacuum 
cleaner is suitable, and necessary, for this) to near zero 
superpressure in order to note whether the volume change 
persisted. It was found that the tetroons are semi-elastic 
with about one-half the volume change remaining after 
deflation to a small superpressure. Figure 2 shows, by 
small circles, the 42-in. tetroon volume change as a 
function of superpressure when the tetroons are first 
brought up to a superpressure of 150 mb. The analytic 
expression approximately fitting these circles is 


AV =0.0038(e° °?— 1) (5) 


where AV is the volume change in cubic meters and Ap 
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Figure 2.—Volume changes of 42-in. tetroons as a function of 
. . . 
tetroon superpressure (circles) and the fit ofgan_analytie expres- 
sion to these values (dotted line). 
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line) and the volume-change term for various pressure heights 
(dashed line) in equation (3) as a function of tetroon super- 
pressure (abscissa). The ordinate gives the change in the weight 
to be attached to the tetroon to enable it to fly at a different 
pressure surface Ap millibars away. 


is the superpressure in millibars. Of great interest are 
the conditions under which the volume change of the 
42-in. tetroon given by equation (5) is of sufficient mag- 
nitude to over-compensate the density-change term in 
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TABLE 1.—Tetroon tracking experiments during 1958-59 





Location | Date Tracked by | Max. range Results 
| (n. mi.) 

Naval Research Labora- | 958 | Mark 25 (3-em.), 19 Good 

tory, Chesapeake Bay ANA/SPQ-2X1 

Annex. (10-cm.) fire control | 

radar and optical 
theodolite. 

Oak Ridge, Tenn. (Oak 9-58 | Decca 40 (3-cm.) |? Ambiguous, 


Ridge National Labo- 
ratory). 

Dawsonville. Ga. (Geor- 
gia Nuclear Labora- 
tory, Lockheed Air- | 
craft Corp.). | 


Weather radar 
(fixed antenna). 
10-58 | Decca 40 (3-cm.) ? 
| Weather radar 
(fixed antenna). 


Ambiguous, 


} 
| | 
Idaho Falls, Idaho (Na- | 11-58 | Modified APS-3 (10- | 4 Poor. 
tional Reactor Testing | em.) radar. 
Station). | | 
Oak Ridge, Tenn_-__-__-- | 3-59 | FPS-10 (10-cm.) radar | 4-6 | Poor, no 
| and optical theodo- | radar 
| _ lite. } | _ return. 
Manassas, Va. (U.S, Air | 3-59 | FPS-3 (23-cm.) radar _| (Altitude Fair. 
Force, Air Defense 5,000 ft.) 15, 
Command). (Altitude | Good. 
| 10,000 ft.) | 
| 60. | 
Hatteras, N.C. (U.S. | 959 | Modified SP-1M (10- | 92 --| Excellent. 
Weather Bureau sta- | to em.) radar, | | 
tion). 


10-59 | 





equation (3), since if this criterion is satisfied the tetroon 
is no longer a constant level balloon. Figure 3 shows 
the relative magnitudes of the density-change term (solid 
line) and the volume-change term for various pressure- 
heights (dashed lines) in equation (3) as a function of 
tetroon superpressure (abscissa). Since the volume- 
change term in figure 3 is plotted as the negative of the 
way it appears in equation (3), it is seen that if the tetroon 
superpressure much exceeds 100 mb., the ordinate, 5.Y/ép, 
becomes negative and the tetroon will ascend until it 
bursts. Consequently, if these tetroons are to be flown 
at ambient pressures of less than 900 mb., it is necessary 
to launch them in a partially deflated state. Procedures 
for determining the weight to be added to the tetroon in 
order to attain flight at a given density surface are given 
in the appendix for cases when (a) the tetroon volume is 
assumed constant, (b) the flight is to be made at low 
levels and the tetroon volume is assumed to vary according 
to equation (5), and (c) the tetroon is launched in a semi- 
deflated state for flight at high level and the tetroon 
volume is assumed to vary according to equation (5). 


4. SUMMARY OF TETROON FLIGHTS DURING 1958-59 


Table 1 gives a summary of tetroon flights made 
through the year 1959. Tetroon flights were first made 
at the Chesapeake Bay Annex of the Naval Research 
Laboratory, where the attachment of a metalized mesh to 
the aluminized tetroon permitted the balloon to be 
tracked for 9 mi. Without this mesh the radar returns 
were sporadic, probably depending upon whether the face 
or the corner of the tetroon was directed toward the radar. 
Figure 4 shows three altitude-versus-time plots of the most 
successful tetroon flight at the Chesapeake Bay Annex, as 
obtained from radar information and a combination of 
radar and visual theodolite information. Of interest in 
this figure is the evidence that the tetroon oversho! its 
floating level and subsequently performed damped ver' ical 
oscillations of a period (in this case) of about 20 min 








us, 


us. 


on 
Ws 
lid 
ire- 
_ of 
me- 
the 
oon 
dp, 
| it 
wn 
ary 
ures 
1 in 
ven 
e is 
low 
ling 
‘mi- 
‘oon 


-59 


lade 
iade 
arch 
h to 

be 
urns 
face 
dar. 
nost 
x, as 
n of 
t in 

its 
‘ical 





i960 


JULY 


Laier tetroon flights at Oak Ridge, Tenn. and Dawson- 
ville. Ga. were not successful. Conditions for tracking 
were poor with precipitation echoes scattered over the 
Decca PPI scope. In addition, the ground clutter was 
bothersome at both stations. Poor results also were 
obtained from the tetroon flights made at Idaho Falls, 
Idaho, in November 1958, probably owing to the low 
power of the modified APS—3 radar. 

The Air Defense Command (ADC) radar at Oak Ridge 
was unsuccessful in attempting to track tetroons in March 
1959. Again precipitation was present in the vicinity of 
the radar and it is believed that at least one tetroon 
entered a region of precipitation and was forced to the 
ground. On the other hand, tetroon flights made from 
the U.S. Air Force ADC site at Manassas, Va. during this 
same month were moderately successful, and sporadic 
tracking of a 10,000-ft. flight was carried out for a distance 
of 60 mi. A longer track might have been obtained except 
for the fact that the constant level tetroons tend to fly out 
of the cone of vision of the upward-directed ADC radars. 
The long strips of aluminized Mylar suspended from the 
tetroon corners on some of the Manassas flights did not 
noticeably improve the radar reflectivity of the balloon 
system and proved most unwieldly during launching. 

Finally, in September and October 1959, tetroon flights 
were made from the Weather Bureau Station at Cape 
Hatteras, N.C. Tracking was accomplished by means of 
the manually operated SP-1M radar of the Weather 
Bureau. These flights were most successful, and conse- 
quently, the remainder of this article is devoted to these 
flights and the meteorological results obtained therefrom. 


5. TETROON FLIGHTS FROM CAPE HATTERAS 
A. EQUIPMENT AND PROCEDURES 


The SP-1M radar at Cape Hatteras has a power of 750 
kw., a wavelength of 10 cm., and a conical beam width of 
3.4°. The range accuracy of the SP-1M is stated only 
as being within 200 yd. If it is assumed that this means 
that 99 percent of the time the error in range is less than 
200 yd., then assuming a Gaussian distribution of range 
errors, the average range error would be 77 yd. Since the 
error in the difference between two quantities equals the 
square root of the sum of the squares of the individual 
errors, the average error in distance between two radar- 
determined ranges would be 109 yd. Thus it is estimated 
that the average SP-1M-determined tetroon speed error, 
AV, in knots, for a tetroon moving radially away from 
the radar would be given by 


AV=3.2/T (6) 
where 7’ is the time interval between positions in minutes. 
One inight like to make 7 large in order that the error in 
letroo’) speed determination be small. On the other 
hand the larger 7’ is chosen, the more the high frequency 
speec )scillations are damped. The compromise adopted 
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Ficure 4.—Tetroon altitude as a function of time 
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Flight 2 from the Chesapeake Bay Annex. 
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Figure 5.—Theoretical average errors in tetroon speed (solid line) 
and height (dashed line) as a function of range for the SP-1M 
radar at Hatteras. The remaining pair of lines gives the observed 
average absolute deviation of speed and height from 20-min. 
average values as a function of range for the tetroon flights 
from Cape Hatteras. 


here involved the determination of speed over a 5-min. 
interval, thus presumably yielding an average speed error 
of 0.64 kt. and introducing a 50 percent reduction in the 
magnitude of oscillations of 8-min. period. Figure 5 
shows by means of solid lines this theoretical error in wind 
speed and, for comparison, the regression line of the ob- 
served average absolute deviation of tetroon speed from 
the 20-min. average value as a function of radar range. 
Since there is no reason to expect the actual speed oscilla- 
tions to increase in magnitude with distance from the 
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radar, the deviation of the observed line from the hori- 
zontal probably indicates that actually there is a greater 
uncertainty in the wind speed determinations as the radar 
range increases. 

Large errors in azimuth and elevation angles may result 
from the use of the SP-1M radar because of the 3.4° coni- 
cal beam width (defined as the angular measure of the 
distance from the power P at the center of beam to 0.5 P) 
and the lack of a precision target indicator (PT). 

If, in the absence of the PTI, it is assumed that it is 
equally probable the tetroon will be positioned anywhere 
within the beam the maximum angular error would be 
1.7° with an average error of 0.85°. The average SP-1M- 
determined height (or lateral) error A//, in feet, would be 
given by 

AH=90R (7) 
where ? is the range in nautical miles. Figure 5 shows by 
means of dashed lines this theoretical error in height de- 
termination and, for comparison, the regression line of the 
observed average absolute deviation of tetroon height 
from the 20-min. average value as a function of radar 
range. The observed deviation increases much more 
slowly than the theoretical, presumably because of the 
necessity for aligning the target more nearly in the center 
of the radar beam in order to get a good signal return at 
great ranges and the subjective tendency not to vary the 
elevation angle of the radar at these ranges unless forced 
to do so. The tetroon targets were manually tracked, 
with one man operating elevation and azimuth controls 
and one on the range scope. The balloon position was 
read from these dials at 1-minute intervals. In any event, 
it is obvious that the vertical and cross-stream motions of 
the tetroon obtained from the SP-1M radar should be 
treated with caution. 

One of the purposes of the Cape Hatteras tetroon flights 
was to test the improvement in radar return resulting 
from use of a metalized nylon mesh obtained from Suchy 


0130 Oct. 1 


Surface pressure patterns at 12-hr. intervals during the period of tetroon flights from Cape Hatteras. 
Cape Hatteras Station is indicated by the heavy dot. 


1330 0130 Oct. 2 


The position of the 
Isobars drawn at 4-mb. intervals. 


TABLE 2.—Tetroon flights from Cape Hatteras, N. C. 


' 


| 
Time of 


Flight | Aver- Aver- 
Flight | launch- dura- | Max. age oh age or Flight 
no. ing Date tion | range | height) (100’s | speed | (kt.) | equipment 
(EST) (hr.) |(n.mi.)} (100’s | of ft.) | (kt.) 
of ft.) 
1 1108 | 9-30-59 My 52 ll (*) 30.6 (*) 42-in. tetroon. 
2 1422 | 9-30-59 4.0 92 | 42 22 23.0 3.3 | 42-in. tetroon 
| with 10-em. 
| radar mesh. 
3 2037 | 930-59 1.9 5D 16 ‘ 26.3 2.7 | 60-in. tetroon 
with 3-cm. 
} radar mesh. 
4| 0825 | 10-1-59 | 4.4 60 40 16 13.6 1.7 | Train of 3 42- 
in. tetroons. 
5] 1413 60-in. tetroon. 


10-1-59 | 5.0 |} 51 143 22 10.2 2.3 


*Not computed—tracking sporadic. 
o.=Standard deviation of height. 
o,=Standard deviation of speed. 


Division, Inc. of New York. This mesh was formed into 
a evlinder and hung like a skirt from the mid-section of 
the tetroon. The mesh weighed only 50 gm. and _ pro- 
vided an additional radar-reflecting area of 30 ft.*. Fur- 
thermore, the mesh skirt, hanging in the form of a cylinder, 
greatly lessened the problem of signal fading associated 
with the shape of the tetroon. The success of the mesh 
can be judged from the fact that the two tetroons on 
which it was installed were tracked over the radar horizon. 


B. TETROON TRAJECTORIES 


Figure 6 shows the surface pressure pattern at 12-hr. 
intervals during the period of tetroon launchings from 
Cape Hatteras. The flights were made (inadvertently) 
while hurricane Gracie was in the process of being trans- 
formed into an extratropical storm over the Middle 
Atlantic States. Consequently, at least 5 min. tracking 
time was lost each hour due to use of the radar for synopti¢ 
observations. Table 2 gives information on the five 
successful tetroon flights made during this period. From 
comparison of the tetroon release time, as presented i 
the table, and the synoptic maps one can visualize the 
synoptic situation during the tetroon flight. It is see 
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from table 2 that in addition to tetroons with metalized 
mesh attached, tracking was also performed on a plain 
42-in. tetroon, a plain 60-in. tetroon, and on three 42-in. 
tetroons tied together. 

Figure 7 shows the trajectories of the five flights. The 
tetroon positions are plotted at 10-min. intervals while 
the tetroon velocity, determined from the distance and 
direction between successive positions, is indicated in 
conventional meteorological form (in knots). The num- 
bers beneath the positions indicate the tetroon height to 
the nearest 1,000 ft. All flights except Flight 5 were 
ballasted to float at 900 mb. (about 3,000 ft.), while 
Flight 5 was ballasted to float at 700 mb. (about 10,000 
ft.). The error in altitude determination of Flight 5 
resulted from the assumption that equation (5) applied 
to 60-in. tetroons as well as to 42-in. tetroons for which it 
was derived. It is noted that near the end of the tra- 
jectories on Flights 2 and 4 there was a sudden jump in 
the tetroon height as indicated by the radar. Since it is 
most unlikely that the tetroons appreciably changed 
altitude, it is probable that these height jumps depict the 
complex refraction and reflection phenomena of the radar 
beam near the radar horizon. Therefore, while the hori- 
zontal trajectory of the tetroon may be determined near 
to and even beyond the radar horizon, it is probable that 
the vertical positions of the tetroon indicated by the 
radar near the radar horizon are of little significance. 
If, however, the variations in signal are random, spectral 
analysis can separate significant vertical motions. Even 
if the signal variations are ordered there is the possibility, 
by comparison with data closer to the radar, of identifying 
the spurious variations due to anomalous propagation. 


C. VELOCITY FLUCTUATIONS 


Figure 8 shows the speed as a function of time as ob- 
tained from the four accurately positioned tetroon flights. 
This speed was determined only from the change in radar 
range apropos the discussion in 5a concerning the in- 
accuracy of azimuth angles determined by the SP-1M 
radar. The dashed speed traces indicate places where 
the tetroon speed had to be interpolated, either due to 
use of the radar for synoptic purposes or to loss of con- 
tact with the tetroon through a fading radar signal. 
The vertical arrows near the beginning of each trace 
indicate the approximate time the tetroons reached flight 
altitude. Most apparent in figure 8 is the long term in- 
crease of wind speed with time on Flights 2 and 3. A 
glance at the first 1330 est map in figure 6 shows that this 
speed increase is in agreement with the increase of surface 
pressure gradient as one moves north of Cape Hatteras. 
On Flight 5 there was a most surprising periodicity in 
wind speed near the beginning of the flight. Several 
times the wind speed changed by 6 kt. in about 13 min. 
of time (2.8 n. mi. travel distance) with some evidence 
that ‘he speed increase was more abrupt than the speed 
decrexse. More will be said about this periodicity later. 
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The significance of the speed changes shown in figure 8 
can be estimated from the discussion of the range accuracy 
of the SP-1M radar and the deviation of the wind speed 
from the 20-min. averages presented in section 5a. 


For each of the four well-positioned tetroon flights, the 
contribution of oscillations of various frequency to the 
variance of the series was determined for the wind speed, 
V, the along-stream, V,, and cross-stream, V,, velocity 
components (where the “stream” is defined by the total 
wind vector for a run) and the vertical motion, W, in 
the manner indicated by Tukey [14]. Spectra were de- 
termined separately for both the earlier portion of the 
flight, when little interpolation of the velocity values 
was required, and for the entire flight. 
were then averaged so that the “good” portion of the 
flight was weighted twice as heavily as the ‘poor’ 
portion of the flight. While considering the spectral and 
aucocorrelation curves so obtained it should be kept in 
mind that the “predominant’”’ periodicities which result 
are to some extent related to the interval and length of 
sampling, and furthermore that the fraction of the energy 
of the velocity represented by V, and V, at a given 
frequency is in part a function of the orientation of the 
coordinate axes. The degree to which these factors bias 
the results presented herein will be more readily apparent 
upon analysis of the precise trajectory data obtained by 
means of tetroon flights from the Wallops Island station 
of the National Aeronautics and Space Administration. 


These spectra 


The variation with frequency of the speed variance 
for the individual flights, as well as the mean variance 
for all the flights, is indicated in the top diagram of figure 
9 for periods of oscillation varying from 3 hr. to 6 min. 
The greatest speed variance is associated with the low 
frequency oscillations, which are probably “synoptic’’ or 
inertial in character. Flights 5 and 3 yield a pronounced 
tendency for speed fluctuations of period near 26 min. 
These two flights have sufficient influence on the mean to 
make the next most pronounced mean periodicity one of 
26-min. period. Less pronounced peaks in the mean 
spectra occur at periods of 12-13 min. and 8 min. Since 
the speed has been averaged over a 5-min. intervai, the 
spectral peak at 12-13 min. might be due to errors in 
range determination. 

The variation with frequency of the variance of the 
vertical tetroon motion is given in the lower diagram of 
figure 9. With the exception of minor peaks of dubious 
significance, the mean variance of tetroon vertical motion 
is at a maximum near a period of 13 min. The three 
possibilities with regard to this peak are: (a) the 13-min. 
period in tetroon oscillation represents a similar period 
in vertical air motions, (b) the 13-min. period reflects 
only vertical tetroon oscillations and does not reflect the 
vertical air motions at all, and (c) the 13-min. period does 
not even represent the predominant period of tetroon 
oscillations but is indicative only of errors in the elevation 
angles obtained from the SP-1M radar (see section 5a.). 
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Figure 7.—Trajectories of the five tetroon flights from Cape Hatteras. The tetroon positions are plotted at 10-min. intervals, while the 
tetroon velocity (in knots) is indicated in the conventional meteorological form. The tetroon height in thousands of feet is indicated 
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Wind speed as a function of time after release as 
obtained from the four accurately positioned tetroon flights 
from Cape Hatteras. The dashed traces indicate interpolated 
speed values. The arrows indicate the approximate time the 
tetroons reached flight altitude. 


Figure 8. 
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Figure 10. 
quency of oscillation for individual tetroon flights and the mean 
of all flights (solid line) from Cape Hatteras, 


Cross-stream velocity variance as a function of fre- 


Discussion of these possibilities is reserved for the next 
It is interesting to note, however, that ex- 
cluding the variance associated with oscillations of more 
than a 1-hr. period the flights with relatively large values 
of the speed variance have relatively small values of the 
vertical motion variance. ‘Thus Flights 3 and 5 each have 
a relatively large speed variance (1.9, 4.1 kt.*) and rela- 
tively small vertical velocity variance (2.7, 2.8 kt.*) while 
Flights 2 and 4 each have a relatively small speed variance 
(1.0, 0.7 kt.*) and a relatively large vertical velocity 
variance (10.0, 6.8 kt.*). In considering these values it 
should be remembered that Flight 3 was an evening 
flight, while Flight 5 was a flight at 14,000 ft. In any 
event, this negative correlation between the magnitudes 
of speed and vertical velocity variances suggests that most 


subsection. 


of the speed variance is not due to vertical oscillations 
of the tetroons in regions of vertical wind shear. 
Figure 10 shows the variance of the cross-stream velocity 


Speed variance (top) and vertical velocity variance 


om) as functions of frequency of oscillation for individual 
n flights and the mean of all flights (solid line) from Cape 


component as a function of frequency. 
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component is very large compared to that of the speed. 
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Figure 11.— Mean phase lag between tetroon vertical motion and 
speed as a function of frequency of oscillation for the Cape 
Hatteras flights. The dashed lines drawn at a phase lag of 90° 
indicate where the maximum upward tetroon motion followed 
the maximum tetroon speed by a phase lag of 90° (14 wavelength). 


This is due to uncertainties in azimuth determination by 
the SP-1M radar. Of the two peaks in the cross-stream 
veloctty variance, the one at a period near 10 min. is 
probably due to azimuth errors while the one at a period 
near 17 min. may reflect a real periodicity in the atmos- 
phere during this time. 

The variance of the along-stream velocity component 
was also determined as a function of frequency. How- 
ever, since it so closely mirrors the results presented in 
figure 9 for the speed, it is not reproduced here. It should 
be mentioned that Mantis [13], from analysis of accu- 
rately positioned superpressured balloon flights at 30,000 
ft., found a minimum of velocity variance at periods of 
12-24 min., certainly not in agreement with the results 
found here. 


D. CONSIDERATIONS OF TETROON MOTIONS AS INDICATORS 
OF AIR MOTIONS 


In estimating the degree to which oscillations in the 
tetroon velocity represent oscillations in air velocity, it is 
desirable to have knowledge of the phase lag between the 
oscillations of various velocity components. This phase 
lag as a function of frequency of oscillation may be ob- 
tained from evaluations of the co-variance and quadrature 
variance in the manner indicated by Van der Hoven and 
Panofsky [16]. 

Let us first consider the significance of the speed 
periodicities indicated by the upper diagram of figure 9. 
It is apparent that if the tetroon does not follow the 
vertical air motion but instead, perhaps due to overshoot- 
ing its flight level initially (fig. 4), undergoes vertical 
oscillations independent of vertical air motions in a 
region of vertical wind shear, then a fictitious periodicity 
in wind speed will be introduced. We should therefore 
be suspicious of any wind speed periodicity which is 
associated with a periodicity in vertical tetroon motion. 
From figure 9 this association is most obvious for oscil- 
lations of 12-13-min. period. For further evidence we 
turn to the phase lag between wind speed and vertical 
tetroon motion (fig. 11). Since, for all four of the tetroon 
flights under consideration, the wind speed was decreasing 
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with elevation, a tetroon performing vertical oscillations 
independent of vertical air motions would show the 
maximum speed preceding the maximum upward motion 
by a phase lag of about 90°. In figure 11 the dashed 
horizontal line indicates this phase difference and it js 
seen that oscillations of 40-min. period and oscillations 
of period 13-8 min. possess this 90° phase difference. 
Therefore the 12—13-min. periodicities in the wind speed 
may well be due to vertical oscillations of the tetroon in 
a region of wind shear. On the other hand, the 26-min, 
periodicity in the wind speed is certainly not due to 
vertical oscillations of the tetroon, since the maximum 
wind speed follows rather than precedes the maximum 
upward motion with a phase difference of about 90° 
for oscillations of 26-min. period. 

The above discussion, while suggesting that variations 
in tetroon speed of 12—13-min. period are due to vertical 
oscillations of the tetroons with this period, does not 
prove that the vertical oscillations were tetroon oscil- 
lations only and not those of the air, since a tetroon 
enclosed in a blob of air would also experience, although 
to a lesser extent, the component phase lags mentioned 
above. Therefore, we next determine whether the pre- 
dominant period of vertical tetroon oscillation is reason- 
ably close to the period which would be expected of an 
air parcel. By relating the vertical acceleration of an 
air parcel to the bouyancy force acting upon it, it can 
be shown that the period of vertical oscillation (r) of 
an air parcel in the atmosphere is approximately given by 


Tr am 
=2r_/——*— 8 
. Wass \ ) 


where 7) is the absolute temperature in the vicinity of the 
air parcel, g is the acceleration of gravity, y, is the process 
lapse rate (usually assumed dry adiabatic) and y is the 
lapse rate [9]. This equation is not exact, since the effects 
of the surrounding air are not considered. These effects 
tend to make the period of oscillation greater, the more so 
the greater is the horizontal dimension of the cell in rela- 
tion to its vertical extent. For example, in cells with the 
dimensions of Bénard cells (cell height three times cell 
diameter), the period of vertical oscillation would be 
1.46 times the period given by equation (8) [5]. Since 
it is impractical to estimate the cell dimensions from 
the tetroon flights it must suffice to say, that if the tetroon 
follows the 3-dimensional path of an air parcel, the period 
of vertical oscillation of the tetroon should be between 
rand 1.57. Table 3 gives, for the 4 flights, the period 
of maximum variance of tetroon oscillation in the vertical, 
as obtained from the lower diagram of figure 9, and the 
period computed from equation (8) by utilizing radiosonde 
ascents at Cape Hatteras to yield information on the mea 
lapse rates and temperatures at the average floating level 
of the tetroons. In all cases the predominant periodicity 
in tetroon oscillation in the vertical exceeded that given by 
equation (8), with the multiplication factor being !.1 02 
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FigurRE 12. 


Positions of maximum tetroon speed (horizontal arrows) and maximum tetroon upward motion (vertical arrows) relative to 


tetroon transversal velocity components (schematic wave-shaped trajectories) for oscillations of various periods (in minutes) for the 
individual numbered flights and the mean of all flights (long arrows) from Cape Hatteras. 


Flights 2 and 5, 1.3 on Flight 3, and 1.8 on Flight 4. 
On the other hand, the natural period of oscillation of an 
absolutely constant volume tetroon (neglecting aero- 
dynamic drag which, for the tetrahedrons, must be a com- 
plicated (and unknown) non-linear function of the restor- 
ing buoyancy force) in an atmosphere with the given 
stability would be about one-fourth the observed 13-min. 
period of vertical oscillation, while the period of vertical 
oscillation of a completely elastic helium-filled balloon 
(helium lapse rate equals 1.3° C. per 100 m.) would be 
about one-half the observed predominant period. Thus 
the predominant period of 13 min. is not due to natural 


TABLE 3.—Comparison of pre- 
dominant period of tetroon ver- 
tical motion and _ theoretical 
period of vertical air motion for 
Cape Hatteras flights 


Predominant | Period computed 
Flight | tetroon period from equation 


(min.) (8) (min.) 
2 10.6 9.4 
3 15.0 | 11.2 
4 13.8 | 7.7 
5 | 8.6 8.0 
| 
| 





vertical oscillations of the tetroon and the hypothesis 
that the tetroon at least partially follows the vertical 
air motion is certainly not contradicted except perhaps 
in the case of Flight 4. However, since this latter flight 
consisted of 3 tetroons tied together (which probably had 
different natural floating levels), it would not be surprising 
to find vertical air motions poorly represented by the 
vertical motions of this tetroon train. 


E. PHASE LAG BETWEEN CROSS-STREAM FLOW, 
AND MAXIMUM SPEED AND UPWARD MOTION 
FOR OSCILLATIONS OF DIFFERENT PERIODS 


Figure 12 gives the position of maximum wind speed 
(horizontal arrows) and maximum upward motion (ver- 
tical arrows) along schematic wave-shaped trajectories 
for each flight, and the mean of all flights, for various 
periods of oscillation. In the mean for all flights it is 
seen that for most oscillation periods the speed was at a 
maximum ahead of the trough in the trajectory, the only 
exception being oscillations of a 45-min. period. Thus, 
if it is assumed that the large-scale air flow is nearly 
geostrophic and that small-scale geostrophic oscillations do 
not exist, then the tetroon moved faster when it was 
moving down the pressure gradient than when it was mov- 








246 


TaBLe 4.—Autocorrelation coefficients of tetroon speed 
(V), vertical motion (W), cross-stream (V,) and along- 


stream (V,), velocily components for Cape Hatteras 
flights 

; : : : |} 1—R(V,) 
Minutes R(V) RCW) | RCV») R(V,) - 
1— FAC.) 

1 0. 90 0. 61 0. 69 0. 87 0. 42 

2 . 51 | . 32 . 40 . 76 40 

3 72 12 | 15 . 67 . 39 

4 | 62 —.10 —.10 . 57 . 39 

5 | 55 —, 20 | —.22 49 42 

6 52 —. 23 | —.27 | 45 43 

7 51 —.14 —.23 45 45 

8 SO 10 —.18 | . 45 | 47 

9 48 —. 08 —.13 - 46 | . 48 

10 | 47 —.07 —. 06 .47 a 

12 | 44 | —.02 01 . 48 | . 53 

14 | 42 02 05 | 45 | 5S 

16 38 . 04 09 | . 39 | . 67 

18 . 33 | —.03 09 | 31 | . 76 

20 . 32 -—. 04 27 . 76 

24 36 —.01 —.10 3 63 

28 . 23 04 01 .19 .82 


ing up the pressure gradient. Thus, there is some evidence 
that kinetic energy (and momentum) is transported down 
the pressure gradient by these small-scale oscillations. 
Moreover, in the mean for all flights, the maximum 
upward motion of the tetroon took place near the trough 
line for oscillations of a period exceeding 18 min. and near 
the trajectory crest for oscillations of a smaller period. 
Thus, looking downstream the long-period tetroon oscil- 
lations were counterclockwise in a plane normal to the 
trajectory, while the shorter-period oscillations were 
clockwise. Insofar as the isotherms are parallel to the 
contours at this elevation, with cold air to the left of the 
flow looking downstream, this means that the long-period 
tetroon oscillations are direct (warm air rising, cold air 
sinking) while the short-period oscillations are indirect 
(cold air rising, warm air sinking). However, 
according to figure 11, the short-period tetroon oscillations 
in the vertical may not be representative of air oscillations 
in the vertical, the evidence for indirect circulations 
should be treated with caution. Further, until data on 
the temperature field are available on a comparable scale 
with the air motions here discussed, comparisons derived 
from synoptic-scale experience must be quite tentative. 


since, 


F. SOME COMPARISONS WITH OTHER LAGRANGIAN DATA 


Mean autocorrelation coefficients of the speed, V, 
vertical velocity, W, and cross-stream, V,, and along- 
stream, V,, components of tetroon velocity for the four 
Cape Hatteras flights are presented in table 4 for time 
lags from 1 to 28 min. Also presented in this table is the 
ratio of one minus the along-stream and one minus the 
cross-stream autocorrelation coefficients. It can be shown 
[10] that for Eulerian space correlation functions in the 
inertial range this ratio would be expected to be 0.6 for 
two-dimensional isotropic turbulence. It is of interest 
that in table 4 this ratio, while fluctuating between 0.4 
and 0.8, does have a mean near 0.6, suggesting a certain 
similarity between Eulerian space and Lagrangian auto- 
correlation coefficients. 
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Figure 13.—Comparison of autocorrelation coefficients of the 
cross-stream velocity component obtained from Cape Hatteras 
tetroon flights, 300-mb. transosonde flights, and the soap bubble 
data of Edinger. (Note that the data are plotted on log-log 
paper, that the ordinate is in terms of one minus the autocor- 
relation coefficient, and that the abscissa is minutes, hours, and 
seconds, respectively, for the tetroon, transosonde, and soap 
bubble data.) 


In figure 13 the value of one minus the cross-stream 
autocorrelation coefficient for the tetroon flights is plotted 
as a function of time on log-log paper. For purposes of 
comparison with other Lagrangian time scales, figure 13 
also shows similar curves for the cross-stream velocity 
components as obtained from transosonde data [3] and 
the soap bubble data of Edinger [4]. Note that the 
abcissa in this figure represents seconds for the soap 
bubble data, minutes for the tetroon data, and hours for 
the transosonde data. It is seen that initially all three 
curves are approximately straight lines with slopes near 
45° on the log-log plot, showing that for all three scales of 
motion 1—R(V,)ot as would be anticipated from La- 
grangian turbulence theory [11]. However, since pro- 
nounced periodicities are present in the cross-stream wind 
component along both the tetroon (17-min. period) and 
transosonde (46-hr. period) trajectories, the above pro- 
portionality between time and 1—R(V,) soon breaks 
down. The soap bubbles could not be tracked for a 
length of time sufficient to delineate the predominant 
periodicity in that small-scale flow. Replacing the sign 
of proportionality by a sign of equality, we find that 
initially the autocorrelation coefficient for the tetroon 
data is given by 


R(V,) =1—4.9t/po (9) 
while for the transosonde data it is given by 


R(V,) =1—3.7t/po (10) 
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where po represents the predominant periodicities of 
oscillation (17 min., 46 hr.) cited above. These values for 
the constants differ from the theoretical Eulerian-time 
constants deduced by Ogura (1.67) and Gifford (1.18), 
and reported by Gifford [6], as would be expected. An 
interesting point is the ratio of these constants, 4.15 for 
the tetroons and 3.14 for the transosonde flights. Hay 
and Pasquill [8] have related Eulerian and Lagrangian 
turbulence statistics (on the assumption that the cor- 
relations decay with time in the same manner) by the 
factor 8. They tabulated 8 values for a large range of 
scales together with the associated turbulence intensity. 
The scale of our experiments is intermediate between the 
diffusion experiments and the air trajectories they ex- 
amined, and computed turbulence intensities for the 
tetroons (table 5) are also of the right order to fit this 
intermediate scale. Thus, these data do not contradict 
the use of a seale factor to relate Lagrangian and Eulerian 
data and indeed support the selection of a 8 value near 4. 

It is also interesting to reverse the usual procedure and 
to work from these Lagrangian statistics to an estimate of 
Eulerian statistics. Following Gifford’s [7] equation (6), 
one can estimate the position of the Eulerian spectral 
maximum. Such estimates from the tetroon flights 
indicate Eulerian maxima in the cross-wind component, V,, 
between 12 and 43 cycles per hour, a frequency not in 
disagreement with the high frequency peak computed by 
Van der Hoven [15] for the Brookhaven tower data. The 
comparison suffers from a great difference in the height 
of the observations but it is, at least, not contradictory. 
The peak in the vertical spectra can be shifted in the same 
manner, giving Eulerian peaks at frequencies from 37 to 
68 cycles per hour. Comparison of these values with those 
obtained by Jones [12] from wind inclination data at 2,000 
ft. again shows reasonable agreement. (However, Jones’ 
data do not extend far enough into the low frequency 
range to completely position his spectral peak.) 

Taking into account the difficulties in positioning the 
tetroons previously discussed, one can examine the total 
variances in relation to any assumption of isotropy. From 
table 5 it is seen that in no case were the variances in the 
three components equally distributed and the lateral, V,, 
to vertical, W’, ratio ranged from about 2 to 4. Flight 3 is 
particularly interesting in this regard. It was the only 
nocturnal flight and indicates both reduced turbulence 
and a smaller V,/W ratio. 


6. CONCLUSION 


The foregoing shows that data of interest and im- 
portance can be obtained from low-level tetroon flights 
even when the radar tracking is not too precise. During 
January 1960, eight tetroon flights were tracked by the 
excellent FPS-16 radar at the Wallops Island station of 
the National Aeronautics and Space Administration. 
The  PS—16 radar positions the tetroons with a root mean 
square error of only 5 yd., and for these flights radar 
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TABLE 5.—Variance distribution and turbulence intensilies 





| | 





| Total variance Average Turbulence intensity 
(kt.2) wind - Z * ste 
Flight No. Pie - speed | 
(kt.) v7 [Tv | Tie, ce —_—/) 
2 2 2 " | 2 Ow!iV | View 
o g,, * c.,” _ | 
iP ' A u , 
ae 40.8 4.0 | 10.6 23.0 0. 28 | 3.6) 0.14 7.1 
3.. 6.9} 7.8) 3.1 26.3 10 | 10.1 | 07 14.9 
5. 18.7 5.6 4.9 10, 2 42 | 2.4 22 4.6 
Average -| 25.2) 5.4 6.9 "17.7 . 28 3.5 16 6.8 





*Weighted Average. 


range, azimuth, and elevation data were available at 
10- or 30-sec. intervals. Consequently, information can 
be obtained from these flights on high frequency Lagran- 
gian oscillations that could not possibly be obtained from 
the SP-1M radar. Thus, with the inclusion of transosonde 
data, information becomes available (admittedly at 
different elevations) concerning Lagrangian fluctuations 
with periods from a few tens of seconds to a week, with 
only a slight gap at periods of 3-8 hr. It is hoped that by 
passing tetroons from one radar station to another, 
reliable information can also be obtained on the Lagrangian 
fluctuations in this range. Such flights are next on the 
agenda. In addition, work is progressing on the develop- 
ment of a very light-weight transponder which can be 
attached to the tetroons in order to increase the tracking 
range. The future of the tetroon in air pollution studies, 
and other studies where knowledge of the Lagrangian wind 
fluctuations is essential, appears bright. 
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APPENDIX 
A. Weigh-off procedure assuming tetroon of constant 3. Determine what volume of helium (V,) would 
volume fully inflate tetroon at flight level plus Ap, where 
1. Determine surface tetroon volume (V,) from Ap is less than 100 mb.: 
equation (4). ’ 
2. Determine density surface (p,) at which flight is to V,—VsPatse 
be made. Pas 
3. Determine helium density (p,,) at earth’s surface ; : . = 
; : : ‘ where p,; is surface air density and pays, 1s alr 
from equation of state for helium. ‘Sib got CPE 
4. Weight to be added to tetroon equals V,(p.— prs) — Genaity when — fully aenabes. 
. ye ; 4. Deflate tetroon until free lift of tetroon equals 
W, where W, is weight of tetroon. 
B. Weigh-off procedure for low-level tetroon flights Vile in be, 
assuming tetroon volume varies according to equation ; 
(5) where p,, is surface helium density and W, is 
1. Determine surface tetroon volume (V,) from tetroon weight. 
equation (4). 5. Determine tetroon volume at flight level (V_) from 
2. Determine density surface (p,) at which flight is 
to be made. V.= V,+ 0.0038 (e° °'54? — 1) 
3. Determine tetroon superpressure (Ap) at density 
surface of flight assuming tetroon released in 6. Weight to be added to tetroon equals Vp. 
fully inflated state. Viens— W 
4. Determine tetroon volume at flight level (V_) from Notes: 


V.=V,+ .0038(e* 54? — 1) 

Weight to be added to tetroon equals V.p,— Vspns— 
W, where p,; is surface helium density and W, 
is tetroon weight. 

(. Weigh-off procedure for high-level tetroon flights 
(tetroon partially deflated at the earth’s surface) 
assuming tetroon volume varies according to equation 
(5) 

1. Determine = surface 
equation (4). 

2. Determine density surface (p,) at which flight is 
to be made. 


~T 


tetroon volume (V,) from 


1. The selection of a desired flight altitude cannot be 
entirely arbitrary for the ‘“‘deflated’’ tetroon 
(Method C) since enough free lift (contained 
helium) must be retained to escape surface 
turbulence and clear obstacles in the launching 
area. 

2. There is evidence that marked temperature inver- 
sions can inhibit, or even prevent, a “deflated” 
tetroon from achieving the chosen flight level by 
radically altering the air vs. helium density 
relationships compared to surface inflation con- 
ditions. This phenomenon is under study. 
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AN EXTENSION OF A TABLE OF ABSORPTION FOR ELSASSER BANDS* 


D. Q. WARK AND M. WOLK 


U.S. Weather Bureau, Washington, D.C. 
[Manuscript received July 19, 1960; revised July 25, 1960] 


Elsusser [1] formulated the absorption by a band con- 
sisting of Lorentz-shaped lines equally spaced and of equal 
strengths. 
interval D, the spacing between the lines, and through a 
homogeneous gas of path length wu, to be 


*y 
A=sinh 3| exp(—y cosh 8) Jo(iy)dy (1) 
0 


where 8=27a/D 
and y= (Su)/(D sinh 8) 


in which @ is the half width of the lines at half maximum 
and S is the strength (integrated absorption coefficient) 
of an individual line. 
Kaplan [2] has shown that equation (1) may be repre- 
sented by 
A=sinh 6 exp(—y) [ Jovin) any" + i Min Sha | (2) 
n=1 n=1 
in which a,;=1 
b,=—6/sinh 8, cosh B<3 
nb ns 1 =b,+4, 
Nn =An_-1+6,-1+(—1)"*! c*'/(n—1)! 
c=cosh B—1 


where 


From this expansion Kaplan prepared a table of frac- 
tional absorption, in which 8 ranges from 0.01 to 1.0 and 
y from 0 to 40. 

Radiative transfer problems arising in meteorology 


*This work was supported by the National Aeronautics and Space Administration. 


He found the mean absorption over a spectral ° 


require an extension of the table. In the present work, 
therefore, the series solution of Kaplan has been employed 
to carry the table to the limits normally encountered in 
the atmosphere. Calculations were performed for 8=1.0, 
0.1, 0.01, and 0.001. 
were required in the evaluation of the coefficients. Ex- 
ponential functions of the desired accuracy were obtained 
from Van Orstrand’s [3] tables. 

The ranges of y were 0 to 150,000 or to the point where 
A became unity to four decimal places. For intermediate 
values of 8, an interpolation scheme was devised to 
determine the difference between the correct values of A 
and the error function approximation, 


As many as 60 decimal places 


A=erf (Byy/2 ) (3) 


Table 1 shows the results of this work. The calculated 
values of A are exact to the accuracy shown. The 
interpolated values should be correct to within 0.0001 in 
most cases. All values of A lying outside this table, 
except for 8>1.0, can be calculated precisely to at least 
four decimal places from the error function approximation. 
Kaplan’s table has been included as a part of table 1. 
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ABSTRACT 


Winter-month data obtained at United States Antarctic stations during the IGY are used to compute coefficients 


for regression equations relating the 100-mb. temperatures and heights to the 50-mb. heights. 


It is found upon trial 


that the layer thickness can be estimated as accurately from the 100-mb. temperature alone as from both the tempera- 
ture and 100-mb. height; therefore, the simpler equations in only one variable are used. The average errors in 
computed heights are less than 2 dekameters and the extreme errors 6 dkm. or less when tests are performed on inde- 


pendent data. 


A measure of the average lapse rates within the layer between 100 mb. and 50 mb., derived from the same data as 
the equations mentioned above, proves to be as good a predictor as any of the regression equations. These mean 
lapse rates show a regular variation with both latitude and month. The largest values (greatest instability) occur in 
the coldest months and the highest latitudes, and the greatest layer stability in lower latitudes and warmer months. 

Additional 50-mb. heights are extrapolated from the observed 100-mb. data for all of the stations considered for 
radiosonde ascents reaching 100 mb. but not 50 mb.; means of the computed heights are compared with the means of 
the observed 50-mb. heights to determine if any bias in the observations can be demonstrated. 


1. INTRODUCTION 


During the International Geophysical Year in the 
Antarctic, regular radiosonde ascents at 12-hour intervals 
were made at, inter alia, the United States stations at the 
South Pole, Byrd, Little America, McMurdo Sound, 
Ellsworth, and Wilkes, and at the joint New Zealand- 
United States station at Hallett. In the colder six 
months of the year, April to September, and in particular 
during the four winter months from May to August, a 
large percentage of those ascents were terminated between 
100 mb. and 50 mb., usually due to balloon bursts. The 
50-mb. height data can be considerably augmented by an 
extrapolative technique to obtain reasonably accurate 
values of 50-mb. heights from the 100-mb. data. 

The purpose of this investigation is to find the most 
feasible method of obtaining additional 50-mb. data by 
extrapolation from observed 100-mb. data, and to ex- 
amine the thermal structure of the 100 mb. to 50 mb. 
layer and its variations with latitude and month over the 
Antarctie during the period considered. 


2. EXTRAPOLATIVE METHODS 


The first approach to the problem was to find coefficients 
for the regression equations relating the 50-mb. height, 
As, to the 100-mb. height, Ajo, and the 100-mb. tempera- 
ture, To. This is the same procedure followed by pre- 
Vious investigators in the Northern Hemisphere, as for 


_—_— 





*p 


presented at the Antarctic Symposium of Buenos Aires, November 1959. 


example, Hering und Antanaitis [1], and currently in use 
by the U.S. Weather Bureau Stratospheric Analysis 
Project. 
100-mb. data, the 50-mb. height is completely determined 
by the known 100-mb. height and the thickness, /7,, of 
the 100 mb. to 50 mb. layer. 
equations, 


Since in all cases we are to work from known 


The original regression 


Hyp= A\ Ho ArT t Ay (1) 


may then be transformed to equations of the form 


Ha= Aj H+ A2T i+ Ao, (2) 
where Aj=A,—1. 

The coefficients for equation (2) were computed for 
Little America for the six months April to September 
1957 (228 ascents) to give a range through the period, 
and for the month of June 1957 for Little America, 
South Pole, Byrd and Wilkes (131 ascents) to give a 
range in latitude. The resulting A, coefficients showed 
a range of from —0.05 to +-0.015, and none was signifi- 
cantly different from zero. 
regression equations was obtained from this same data 
by computing the A, and A) coefficients for the general 
equation 


Therefore, a new set of 


Hg=AjT T An (; 


~~ 


A test of equations (2) and (3) on a limited amount of 
independent data (26 cases) resulted in an average error 
of 13 gpm. in the 50-mb. height using (2) and 11 gpm. 
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TABLE 1.—Mean values of thickness of 100 mb. to 50 mb. layer (Ha) 
and 100-mb. temperature (Ti) for winter months in the Antarctic. 
Thicknesses are in geopotential meters and temperatures in degrees 
Celsius. 





i. 








a 
Little | 
South Byrd America | Hallett | Wilkes Monthly 
Pole (80° 8.) v— | (72° 8.) | (66° 8.) mean 
(90° 8.) | Ellsworth | 
| (78°Ss.) | 
| 
| | | | | 

April Ha |--- 4291 4324 4384 4432 | 4343 
Ti | —58.6 —57.5| —55.0] —53.1 | —56. 6 

2 ES 1c connection . 
May Ha | 4064 4081 4122 | 4266 4350 4164 
T 00 —69, 1 —68.3 —66.4|) —59.6|] —57.3 | —64. 6 
June Ha 3943 | 3977 4032 4114 | 4275 4050 
T1100 —74.8| —73.1 —71.1 | 67.6 —0.4 —70.2 
July Ha | 3852/ 3888 | 3896 | 4009 | 4038 3931 
Tico | —79%1| 77.4} —77.3| 72.2} —71.3] —75.7 
August Ha 3780 | 3821 | 3874 | 3946 4041 | 3880 
T 00 —83.4 —81.9 —79.3| —75.6 —71.9 | —79.0 
Sept. Ha | 3863 | 3902 | 3961 4113 4156 3988 
Ti —81.0 —79.5 —76.5 —70.3 | —69. 1 —75.7 
Mean (Lat.) Ha | 3969 | 3992 4033 4138 | 4214 (4058. 26) 
T'100 | —74.2 —73.2 —71.4 —66.8 | -—63.9 | (—70.33) 


using (3). On the basis of this comparison, it was de- 
cided to use equations of the form (3) in computing 
regression coefficients for all stations except McMurdo 
Sound, and for all months of 1957 and 1958 (April to 
September) in which at least 8 ascents reached the 50-mb. 
level. The data from Ellsworth Station and Little 
America were combined, since these two stations are 
approximately the same latitude; a detailed listing of the 
data used is given by the parenthetical values in table 4 
in the rows headed “Comp.”. The MeMurdo Sound 
data and unused data from the remaining stations were 
reserved for testing the results; table 3 gives a breakdown 
of the independent data by station and month. 

In the course of the above computations, intermediate 
steps made it convenient to extract from the data the 
mean FZ, and Ty for each station and month (table 1), 
and also to compute some of the correlation coefficients 
between these two parameters. The latter were generally 
quite high, ranging from +0.89 to +0.99 with a median 
value of those computed (18) of +0.95. Since the original 
computations were based on only 1455 cases in all, and 
individual station-months on from 8 to 43 ascents, it was 
thought advisable to test on independent data not only 
for the regression equations but also for equations based 
on the assumption of a constant lapse-rate of temperature 
within the layer. 

If we assume that there is no significant amount of 
moisture within the layer from 100 mb. to 50 mb., then 
the mean virtual temperature and the mean temperature 
of the layer are approximately equal. The general 
equation which relates the thickness, 7, (in geopotential 
meters) and the mean virtual temperature, 7, (in ° K.) 
of a layer of air between two constant pressure surfaces 
Pp; (base) and p, (top) is 


H,=67.442 T, logy (p:/p2) 
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TABLE 2.—Computed lapse rate, °C./km., of 100 mb. to 50) 
mb. layer. 

: ; a - ae 4 ia Re Mer iq a 

Little | 
South Byrd | America| Hallett | Wilkes | Monthly 

Pole | V—Ells- | | mean 
| | worth | 

(90° 8S.) (80° 8S.) (78° 8.) (72° S.) | (66° S.) | 

<—s a 3 : mc Gh EARP HENLE it 
April. ._. -1.5 | —15 —12 | -10 | -08 | -1L4 
May...- —-19 | —19 -1.8 | -—16 —@.7 | —1.65 
en —2.1 —2.1 | —1.7 | -—14 —3.0 | —1.72 
| ee —2.2 —2.0 | —2.0 | -—-1L7 | —1L5 —1.9 
August... .._- —-19 | —1.6 —1.6 -16 | -11 | —159 
September__.._____- —-1.0 | -0.7 | -0.8 —0.1 | +03 —.79 

Mean for station (Lati- | 5 | 
ae ae —1.6 | -1.53; -125 =. 7 —1.4 


as given, for example, in [2]. We may now write for the 
100 mb. to 50 mb. layer, substituting the mean tempera- 
ture, 7;, for the mean virtual temperature, the equation 


H,=20.302 T,. (4) 


If we now make the further assumptions that the lapse 
rate is constant with height within the laver and with time 
over the period of a month at a given station, then these 
lapse rates can be computed from the data given in table 
1. The computed mean lapse rate, (T,), for each station 
and month, as well as the weighted mean lapse rates for 
each month and latitude, are given in table 2. 

For extrapolation of the 50-mb. heights, we may now 
write a set of equations, based on the assumption of 
constant lapse rates, of the form 


lH. 

= 20.302] 273.16+ Tio+—2—* 

“ , | 13.16-F £100T 3990 

This is the same as equation (4) with the term inside the 

brackets representing the mean absolute temperature of 

the layer for a given Ty (in °C.) if the assumption of 
constant lapse rate holds. 


3. TESTS ON INDEPENDENT DATA 


Independent data from the four months, May to August, 
were used to test four different methods of extrapolating 
the 50-mb. heights from the 100-mb. data, as follows: 

Method I used equations of form (3) for a given latitude 
and month applied to data from the same latitude and 
inonth; Method II used equations of the form (3) for a 
given latitude, with the coefficients computed by com- 
bining the data from the four months, May to August, 
applied to data from the same latitude; Method III used 
equations of the form (3) for a given month, with the 
coefficients computed by combining all data for that 
month, applied to all independent data for the given 
month; Method IV was based on equations of the form 
(5), although to facilitate computations a nomogram and 
two graphs based on these equations, as described in the 
next section, were actually used. 

The results of the testing of these various extrapo!ative 
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TasLe 3.—Results of tests of four extrapolative methods on 
independent data. 
— — —— —————— — — = ——————— — — 
| Mean errors (gpm.) in computed 
| | 50-mb. height 
Station Months ‘LS . ee ee ea 
| cases | | | 
Method | Method | Method | Method 
| Bee ae ee, 
= en ee 2 : _ 
75°-90°S. 
ads a Pa IPR Ey | ae | : | 
See May and Aug., 1958___ 15 | 8.8 15.7 18.0 | 10.7 
Byrd. __.| May to Aug., 1958 __.- 97 12.2 12.4 15.1 | 12.2 
Elisworth _..| May to Aug., 1957 A 28 11.8 10.7 13.6 14.0 
McMurdo _---.| May to Aug., 1957 a 165 14.8 12.7 | 14.7 | 12.8 
Mean AS yal eae 14 21 Se 2.5) 14.9] 12.6 
North of 75°S. 
- - : — — 
Hallett._._.---- | May, June and Aug. | 74 | 22.5 | 23.1 25.3 22.6 
| 1958. | ] 
Wilkes | May-Aug. 1958__...-_- 79 19.8 | (Independent for Method I 
only: not included in mean 
| below) 
Mean sgl lal ie are eek aout CS | 15.2 | 


14.6] 16.9] 14.5 


methods of obtaining the 50-mb. height are shown in 
table 3. The best results were obtained by using Method 
IV, although the small difference between this mean error 
and that from Method II hardly can be significant. In 
the area from latitude 77° S. to the Pole, the average 
errors using Method I, II, or IV were quite consistently 
between 10 and 15 gpm., while from Hallett northward 
more limited amounts of independent data indicated an 
average error using the best methods of 20 to 25 gpm. In 
general, the extreme error was about four times the average 
error. 


4, CONSTRUCTION OF NOMOGRAM AND GRAPHS 


When we substitute in equation (5) the values for the 
layer thickness Hg averaged over all stations for all of the 
6 months, 4058.26 gpm., and the overall average Ty, 
—1.44° C./km., it becomes 


H4= 5486.4 + 20.302 Ti, 


where /7, is in geopotential meters and Ty in °C. For a 
given 100-mb. height Ayo, a first approximation to the 
50-mb. height Hj) can be obtained, based on the mean 
overall lapse rate and layer thickness, from 


Hs = Hyp + 20.302 T 00+ 5486.4; (6) 


this equation (6) was used to construct a nomogram, using 
Hy as ordinate and Tyo as abscissa, from which the 
approximate Hs could be read off directly. The curves 
in figure 1 show the correction (gpm.) which must be 
applied to the nomogram value for latitude (fig. 1A) and 
for date (fig. 1B) to get the best estimate of the 50-mb. 
height using Method IV. Since the nomogram itself is 
simple linear one it is not shown here. 


5S. LAPSE-RATE VARIATIONS 


rT . . 
Th. correction curves of figure 1 are also graphs of the 
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Corrections to extrapolated 50-mb. height: (A) for 
latitude, (B) for date. 


Figure 1. 


deviations of the layer lapse-rate from the overall mean 
value of —1.44° C./km., with a correction of 10 gpm. on 
the ordinate scale the equivalent of a difference in lapse 
rate of 0.24° C./km. from the mean. Viewed in this 
light, figure 1A indicates that the stability of the 100 mb. 
to 50 mb. layer, averaged over the 6 months, April to 
September, increases in regular fashion from the Pole 
northward to the Antarctic Circle. Figure 1B shows that 


in general the layer stability averaged over all latitudes 
shows a gradual decrease from April until mid-July, and 
a comparatively rapid 
September. 

Comparison with the figures in table 1 shows that while 
the layer is cooling until mid-July, the stability is decreas- 
ing; this means that during this period the 50-mb. temper- 


increase from then through 
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ature is lowering more rapidly than the 100-mb. temper- 
ature, as already noted by Wexler [3]. The more rapid 
cooling at 50 mb. has been explained in terms of the 
vertical distribution of the atmospheric gases important 
in long-wave absorption, by Moreland [4], whose inves- 
tigation indicates a maximum rate of cooling in this layer 
near 18 km., which is near the 50-mb. winter level. 

The flattening of the curve in figure 1B in early June 
indicates that at this time, the temperature is falling at 
about the same rate at 100 mb. and 50 mb. Any explana- 
tion of this somewhat anomalous behavior is beyond the 
scope of this paper, although it might be noted in passing 
that June temperatures during 1957 and 1958 have been 
anomalous at other atmospheric levels, as demonstrated, 
for example, in Wexler’s [3] study in the section on the 
“kernlose’”’ Antarctic winter surface temperature pattern. 

Finally, it might be well to note that for allof the 6 months, 
April to September, over the whole range of latitude the 
temperature falls with height in the 100 mb. to 50 mb. 
layer, with the exception of the lowest latitude station 
during the month with most sun (Wilkes, September), 
when an inversion began in this layer. 


6. COMPUTATION OF ADDITIONAL 50-MB. DATA 


Table 4 shows the mean observed 50-mb. heights for 
several months and stations, and the mean computed 
height of the 50-mb. level from the radiosonde observations 
terminating between 100 mb. and 50 mb., using Method I 
for the computations. Method I is chosen here because it 
gave the best results of any of the four methods when 
applied to the same stations for which the coefficients 
were derived. It, therefore, seemed likely that Method I 
would give the best results when applied not only to the 
same station but to the very same month which gave 
the coefficients. 

Table 4 also gives the weighted mean of the observed 
and the computed heights for each month and station for 
which computations were performed; i.e., those for which 
data available and which had a minimum of 8 
radiosonde observations reaching 50 mb. Numbers in 
parentheses following the mean heights are the number 
of 50-mb. heights, observed or computed, used in deter- 
mining the mean. The heights are in geopotential meters, 
with 18,000 gpm. subtracted from each mean to save 


were 


space. 

Observers have noted that sounding balloons in the 
Antarctic, as elsewhere, tend to reach greater heights 
when rising through a relatively warm atmosphere than 
when rising through a relatively cold one. <A_ priori, 
therefore, one might expect the mean of the computed 
heights during the Antarctic winter to be lower than that 
of the corresponding observed heights; however, the 
results tabulated in table 4 do not demonstrate conclu- 
sively any such bias, with the computed means lower 
than the observed about two-thirds of the time and higher 
the remaining one-third of the time. 
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Mean observed and computed 50-mb. heights. Heights mm 
geopotential meters minus 18,000. Parenthetical figures are number 
of ascents used in determining means. Station-months with fewer 
than 8 ascents reaching 50 mb. not included 


TABLE 4. 





| 
Y : | : 

Station April May | June July August Septem- 

| er 
te 
Pole Year(s) | 1957 1957, 1958 | 1957, 1958 1957 =|: 1957, 1958 
Obs. 908 (30)) 640 (61) | 388 (34) 308 (10)) 324 (29) 
Comp. 952 (17)| 627 (48) | 264 (68) | 313 (37)| 289 (48) 
Mean | 924 (47)| 634 (109)| 305 (102)! 312 (47)) 302 (77) 

| | | | | 

Byrd Year(s) | 1957 | 1987 | 1957 | 1957 1957 | 1957 
dbs. 1455 (8) 1004 (20)' 859 (29) 448 (24) 416 (27) | 457 (17) 
Comp. 1453 (11)} 1104 (17)) 864 (21) | 372 (25) 440 (18)| 411 (20) 
Mean 1454 (19)| 1050 (37); 861 (50) | 411 (49) 426 (45)| 432 (37) 

Ellsworth Year(s) | 1957, 1958 | 1958 = i ara 1957 
bs. | 1569 (52)) 1088 (34) ..| 660 (16) 
Comp. | 1690 (10) 9s4 (13) --| 634 (31) 

Mean | 1578 (62)| 1023 (47) | 643 (47) 
Little Year(s) | 1957,1958 | 1957, 1958 | 1957, 1958 | 1957, 1958 | 1957, 1958 | 1957, 1958 
America Obs. 1661 (71)| 1206 (56)| 933 (73) 489 (76) | 424 (61)) 483 (79) 
V Comp. 1709 (22)| 1874 (34)| 942 (18) 455 (31) 416 (28)| 431 (25) 
Mean 1672 (93); 1269 (90); 935 (91) | 479 wa 421 (89) 471 (104) 

| 

Hallett Year(s) 1957 | 1957 1957 on 1957 | 1957 
Obs. 1760 (22)| 1587 (21)| 1138 (13) | 656 (32)| 899 (25) 

Comp. | 1780 (5) | 1511 (21); 1106 (18) | 569 (19)| 900 (24) 

Mean 1764 (27); 1549 (42)/ 1119 (31) | 624 (51)| 900 (49) 
Wilkes Year(s) 1957, 1958 | 1957, 1958 | 1957, 1958 | 1957, 1958 | 1957, 1958 | 1957, 1958 
Obs. 1957 (71) 1742 (63)| 1657 (35) 1030 (38) 850 (35)) 958 (57) 

Comp. | 1939 (7) 1704 (8) | 1642 (13) 942 (50) | 797 (47)| 933 (29) 

Mean | 980 (88) 820 (82)) 950 (86) 


1956 (78); 1738 (71)) 1653 (48) 


7. CONCLUSIONS 


Before listing conclusions to be drawn from this investi- 
gation, it is well to emphasize that only 7 stations’ data 
from two winter seasons were available at the time of 
writing, and that the stations were concentrated in only 
about one-half of the continent. The effect of these 
limiting factors in modifying the conclusions listed is 
discussed in more detail in the results itemized below. 

1. Extrapolated 50-mb. heights obtained using the 
regression equations with the 100-mb. temperature only 
as a predictor were as accurate as those obtained from 
equations using both the 100-mb. height and temperature. 

2. Extrapolation of the 50-mb. using the 
assumption of constant lapse rates yields as good results 
as those of the computed regression equations. 

Because of the limited amount of data, conclusions (1) 
and (2) above should be considered as true only for this 
set of data. It is entirely possible that the small superi- 
ority demonstrated by the constant-lapse-rate method 
and the single-predictor regression equations is due to 
chance, and that use of both the 100-mb. height and 
temperature as predictors might, in a large sample, show 
improvement over either of the other two methods. It 
is further possible that addition of different predictors 
might lead to a further refinement in accuracy; with the 
limited data it would in any case be impossible to demon- 
strate that any such improved accuracy was the result 
of the method’s superiority rather than the result of 
chance, unless the new equations gave 50-mb. heig!its 0 
a rather unlikely degree of accuracy. 

3. The stability of the 100 mb. to 50 mb. layer was, 
in general, least when the mean temperature of the aif 
was lowest; exceptions were noted in early June, whet 
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the s\ability remained constant, although the mean layer 
tempcrature decreased, and in August when some relative 
warming at 50 mb. increased the stability while the mean 
temperature was still falling slightly. 

4. In the mean, there was a lapse of temperature in the 
100 mb. to 50 mb. layer over the portion of the Antarctic 
considered from April through August, and in September 
only lower latitudes showed an inversion. 

Strictly speaking, conclusions (3) and (4) are valid 
only for the specific locations and seasons considered. 
However, these conclusions can probably be extended 
safely to the continental area in which these stations are 
concentrated, although extension to the whole Antarctic 
area would depend on the unsafe assumption of approxi- 
mate symmetry of temperatures about the geographic pole. 

5. From extrapolative methods a considerable amount 
of additional 50-mb. height data accurate on the average 
to about 20 gpm. can be obtained for the winter months 
of the TIGY at the Antarctic stations considered. 
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Weather Note 


RAPID PRESSURE VARIATIONS IN ICELAND 
EYSTEINN TRYGGVASON 


Geophysical Section, Vedurstofa Islands, Reykjavik, Iceland 


[Manuscript received June 27, 1960] 


In the ‘Weather Note’, Monthly Weather Review, 
January 1960 an unusual pressure rise at Yakutat is 
described [1]. As the note asked for information about 
similar or larger pressure rises, I want to give a brief 
description of an unusual cyclone with rapid pressure 
variations that struck the east coast of Iceland on January 
25, 1949. This cyclone had a central pressure of about 
940 mb., its radius was approximately 200 km., and it 
moved with a velocity of about 70 km./hr., N. 15° E. 

The most rapid pressure variations occurred at Dala- 
tangi, where the following pressure observations were 
made: 

Sea level pressure 


January 25, 0900 Gur 977.8 mb. 


1200 973.7 
1500 963.0 
1800 952.5 
2100 975.6 
2400 981.7 


According to the barograph trace at the station (fig. 1), 
the lowest 941 mb., occurred at 1715 Gmr. 
The maximum pressure fall in 1 hour was from 956 mb. 
at 1600 Gmr to 942 mb. at 1700 Gaur, or 14 mb. The 
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Fiagure 1.—Photocopy of section of tarograph trace 
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maximum fall in 3 hours was 26 mb. The maximum 
pressure rise in 1 hour was from 941.5 mb. at 1740 ower 
to 960.5 mb. at 1840 Gat, or 19 mb. The maximum 
pressure rise in 3 hours was 33 mb. 





The barograph trace made at Holar, Southeast Iceland, 
was very similar to that at Dalatangi with minimum 
pressure of 947 mb. and maximum pressure rise of 18 mb. 
in 1 hour. 

Pressure rises of similar magnitude as that described 
in [1] have been observed several times at Icelandic 
weather stations. Three examples observed in Reykjavik 
during the last 12 vears are: 


l-hour pres- 3-hour pres- 
sure rise sure rise 


January 22, 1950 12 mb. 20 mb. 
January 5, 1952 z 18 
March 8, 1953 9 19 

REFERENCE 
1. Mac A. Emerson, ‘‘Pressure Rise at Yaku: at, December 18, 
1959,’ Weather Note, Monthly Weather Review, vol. 88, 
No. 1, January 1960, p. 18. 
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THE WEATHER AND CIRCULATION OF JULY 1960 


Persistent Heat in the Pacific Northwest 


RAYMOND A. GREEN 


Extended Forecast Section, U.S. Weather Bureau, Washinaton, D.C. 


1. HIGHLIGHTS 


July 1960 was especially noteworthy for prolonged, 
excessive heat in the Pacific Northwest. The greatest 
departures of average temperature from normal, more 
than 6° F., were centered on the Oregon-Idaho boundary 
(fig. 1). New temperature records at Boise, Idaho in- 
cluded an all-time high of 111° F. on the 19th, 27 days 
with 90° or higher, 11 days with 100° or higher, and the 
highest monthly average in 62 years. New records at 
other stations appear in table !. Sustained high tempera- 
tures, low humidity, and little precipitation brought, in 
the words of LaMois of the U.S. Forest Service [1] “. . 
extremely critical burning conditions throughout the 
Western States during July. More than 4,000 fires oc- 
curred, most of which started from heavy concentrations 
of lightning storms which struck California, Oregon, 
Washington, Montana, and Idaho. . . . The fire fighting 
bill during July alone was more than $15 million.” 

Abnormal warmth over the western half of the United 
States was accompanied by unseasonable coolness over 
much of the eastern half, with greatest departures in a 
band from northern Texas to eastern Pennsylvania (fig. 
1). New low averages for July were also numerous. 
Cleveland, Ohio headed the list (table 1) with an average 
temperature of 67.6° F. and a departure from normal of 
—§.1° F. 





TaBLe 1.—New record average temperatures for July, and their 
departures from normal, observed in 1960 
Monthly | Departure 
average from nor- 
(° F.) mal (° F.) 

High tem- | 

perature 
Boise, Idaho , 80.7 +5.9 
Nome, Alaska = = 54.3 +4. 7 
Salt Lak: City, Utah 81.2 +4.6 
Pendleton, Oreg 78.3 +3 
Idaho Falls, Idaho 71.8 +2. 6 
Miami, Fla 83.7 +1.1 

Low tem- 

perature 
Cleveland, Ohio 67.6 —6.1 
Youngstown, Ohio 67.9 | —3.7 
Pittsbureh, Pa 68. 6 —3.7 
Philadelphia, Pa *73. 2 —3.1 
Wilmington, Del 2°73, 1 —2.8 

*Equ lin 1895. 
Equaled in 1894. 


Tropical storm Brenda moved up the east coast the last 
few days of the month. It caused little wind damage, but 
heavy rains were reported along its path. 


2. MEAN CIRCULATION 


From June to July high-latitude blocking diminished 
over northern Canada and increased over the Siberian 
Peninsula (+260 feet, fig. 2). During July two 5-day 
mean blocking ridges coalesced with the subtropical ridge 
in the Pacific, and amplification of the circulation occurred 
downstream at middle latitudes. Height anomalies in the 
troughs and ridges of figure 2 were generally larger than 
their June counterparts, and a new mean wave was intro- 
duced into the pattern in the eastern Atlantic. Further 
evidence of large amplitudes in July is given by the dis- 
placement of the principal axis of maximum mean 700-mb. 
winds. In figure 3 it can be seen that in mean troughs of 
the western part of the hemisphere southward displace- 
ment of the “jet stream” occurred, while in the ridges the 
displacement was northward except in the central Pacific. 

The slow-down of the temperate westerlies was not as 
great as one expects to observe with amplification of this 


magnitude. A drop of only 0.2 m.p.s. from the preceding 
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Figure 1.—Departure of average temperature from normal (°F.) 
for July 1960. Noteworthy features include abnormal warmth 
in the Northwest and coolness from the Central Plains north- 


eastward. (From [2].) 
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Figure 2.—Mean 700-mb. contours (solid) and height departures from normal (dotted), both in tens of feet, for July 1960. 


Crreater 


than normal amplitude of mean troughs and ridges in temperate latitudes is suggested by departure from normal pattern. 


month took place in the departure from normal of the 
mean monthly 700-mb. zonal index for the western half of 
the Northern Hemisphere. Five-day mean indices of the 
temperate westerlies completed a minor oscillation about 
the normal by mid-July and remained above normal 
thereafter. 

Large negative changes of mean 700-mb. height anomaly 
from June to July, shown in figure 4, indicated strong 
deepening over the eastern Pacific and the British Isles at 
temperate latitudes, while falls over northern Canada and 
rises over the Asian sector of the Arctic Basin attended the 


migration of blocking at high latitudes. Over western 
Canada, rises reflected the growth of the mean ridge which 
strongly influenced the temperature pattern of the United 
States. 


3. TEMPERATURE 


The mean ridge at 700 mb. over the Pacific Northwest 
was both stronger than normal and persistent. Five-day 
mean height anomalies there remained well above normal 
the entire month. Not once did the weekly average tem 
perature fall below normal at reporting stations in !daho 
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Figure 3.—Mean 700-mb. isotachs (in meters per second) for July 
1960. Solid arrows indicate axes of west wind maxima, dashed 
arrows the normal fer the month. Eastward from the central 
Pacific the primary axis was displaced southward from normal 
in mean troughs and northward from normal in ridges. 


and inland sections of Oregon and Washington [2], end the 
same was true in Montana after the week ending July 3. 
Thus the extent and intensity of warm temperature anom- 
ilies were largely products of persistence. 

One locally effective warming factor was the downslope 
flow indicated by the departures from normal of mean 700- 
mb. height in figure 2. The flow was generally north- 
easterly where centers of maximum departure were ob- 
served in valleys of Idaho, Utah, and California which 
lie southwestward from extensive mountain barriers. 
Perhaps of greater importance was the fact that more than 
80 percent of possible sunshine was realized at most sta- 
tions in the anticyclonic environment of the Northwest. 

Over the remainder of the Nation there was more intra- 
monthly variability of the anomalies of both temperature 
and 700-mb. height. Even so, temperatures averaged at 
least slightly less than normal each week in Pennsylvania 
and sections of Ohio and New York adjacent to Lakes Erie 
wid Ontario. Seattered stations westward as far as the 
Central Plains also reported subnormal averages for every 
week of the month [2]. 

Most of July’s area of below-normal temperature was 
ilso in the below category in the mean patterns for spring 
and for June [3]. Exceptions were a semicircular region 
from southeastern Colorado to central Texas and Atiantic 
Coasta! States northward from Virginia. In general the 
cool te: peratures of July were what might be anticipated 
from te mean height anomaly field and the cyclonic 
Horthw sterly flow in figure 2. In addition, the large 


Ficure 4.—Difference between monthly mean 700-mb. height 
anomaly for June and July 1960 (July minus June) in tens of 
feet. Pronounced deepening occurred in the northeastern sec- 
tions of the Pacific and the Atlantic. Large changes at polar 
latitudes accompanied the migration of blocking from northern 
Canada. 


amplitude of the mean circulation supported a southward 
displacement from normal of the tracks of daily Highs and 
Lows (fig. 5). One branch of the paths of both dipped 
abnormally far into the United States, augmenting the 
southward penetration of cool Canadian air. 

Cooling over the Soushwest from the extreme heat of 
June left that section of the country only slightly above 
normal in July. A similar change took place in the North- 
Changes of the temperature anomalies from June 
to July followed rather closely the changes of 700-mb. 
height anomaly in figure 4. 
along a diagonal band oriented southwest to northeast, 
and general warming over the Northwest and the South- 


east. 


Thus general cooling occurred 


east. Among 100 selected stations the largest change 
was only 2 classes, qualitatively indicative of more than 
usual month-to-month persistence. 
tively verified by the fact that 76 of the stations appeared 
in the 0- or I-class change category, while Namias found 
June to July persistence of these classes 72 percent on the 
average from 1942 to 1954 [4]. 


This is quantita- 


4, TROPICAL STORM BRENDA AND PRECIPITATION 


Brenda formed in the eastern Gulf of Mexico in late 
July. The behavior of the mean circulation prior to and 
during the storm’s short life is depicted in the series of 
overlapping 5-day mean 700-mb. maps in figure 6. On 
the chart for July 21-25 (A) a deep mean trough extended 
along the Atlantic coast well into the Tropics from Lab- 
rador. By July 23-27 (B) the trough had sheared, and 
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Figure 5.—Number of (A) anticyclone passages and (B) cyclone 
passages within quadrilaterals of 66,000 square nautical miles 
during July 1960. Primary tracks are indicated by solid arrows. 
Branches of the tracks of both anticyclones and cyclones pushed 
abnormally far southward into United States. Note the scarcity 


of migrating systems west of the Continental Divide. 


the tropical portion retrograded to Florida, while the 
northern portion of a similarly sheared trough advanced 
from Oklahoma to Indiana. By July 26-30 (C) the 
fractured parts of the troughs were joined with the 
Hudson Bay Low to form a deep mean trough across the 
Great Lakes and the eastern Gulf of Mexico. Brenda 
formed in this trough. 

The subsequent path of the storm coincided rather 
closely with the 30-day mean trough line of figure 2. 
This tendency for storms to move along mean troughs 
has been discussed by Klein [5], who listed numerous 
examples from previous years. With such coincidence 
the question of cause and effect often arises, since it may 
be argued that low pressures in the storm largely deter- 
mine the location of the trough. For example, see the 
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Ficure 6.—Five-day mean 700-mb. contours in tens of feet for 
(A) July 21-25, (B) July 23-27, and (C) July 26-30, 1960. 
Tropical portion of mean trough in the western Atlantie (A) 
retrograded to Florida (B) and joined with sheared portion over 
western Indiana to form deep trough through eastern Gulf of 
Mexico (C) where tropical storm Brenda developed. Heavy 
rains occurred on July 23 and 24 near the mean Low over Okla- 


homa (see A). 
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Selected precipitation totals (inches) for July including 


TABLE 2. : ae 
totals during close proximity of Brenda 





| | ‘Total 
Total for | during 
July Brenda 
ampa, Fla_ ‘ } *20. 59 | 16. 36 
me Fla__ *19. 57 | 10. 79 
Jacksonville, Fla__ } *16. 21 1. 57 
Savannah, Ga_._ 2 *15. 70 | 2.89 
Ft. Myers, Fla-- 13. 76 | 4.42 
Charleston, &.C 11.74 | 4.43 
New York, N.Y----- *9. 97 | 4.90 
Bridgeport, Conn. _- 8.13 3. 57 
Richmond, Va 7. 34 2.39 
Wilmington, Del-- 6. 18 2.51 


*New record for July. 


correspondence between Ramage and Ballenzweig in a 
recent issue of the Journal of Meteorology |6). To throw 
some light on this question, a mean map of the 25 July 
days preceding Brenda was constructed (fig. 7). The 
storm path superimposed on this pattern follows the mean 
trough line rather well. Obviously in this case the low 
pressures of the storm altered the mean 30-day trough 
position very little. 

The storm’s effect on the precipitation pattern (fig. 8), 
however, was more substantial. Rainfall totals from 
Brenda comprised a large proportion of the monthly 
amount at a number of stations, as revealed by table 2, 
where both totals are tabulated for comparison. At 
Tampa, Fla., more than three-fourths of the monthly 
accumulation could be attributed to the storm. In 
general, much higher proportions of the monthly amounts 
were caused by the storm near the area of its formation 
than elsewhere along the path where rapid motion oc- 
curred. During Brenda’s passage new records for 24- 
hour rainfall were established at Orlando and Tampa, 
Fla., and at New York City. 

Most of the precipitation in the Great Lakes Region is 
attributable to frontal Lows traversing the southward 
displaced storm track of figure 5B. An additional feature 
of interest associated with the first of these Lows was the 
outbreak of tornadoes in North Dakota on July 2 and in 
Indiana on the 3d. Other tornadoes were reported later 
in the month, mostly after the 20th, and were most 
frequent from the Southern Plains to Iowa. Only those 
in lowa were clearly associated with frontal troughs, 
while the remainder occurred beneath cyclonic upper flow 
in maritime tropical air. 

Rainfall was more than twice the normal from the Texas 
Panhandle to Arkansas, where moisture was in adequate 
supply in the weak flow from the Gulf of Mexico around 
the mean High centered over eastern Texas (fig. 2). Some 
of this precipitation was frontal in nature, but the heavi- 
est, In eastern Oklahoma and western Arkansas, was 
triggered by a trough which formed during a readjustment 
of mean waves. When the wavelength downstream from 
the retrograding Pacific trough became long, new mean 
'roug!is appeared on the 5-day mean map for July 19-23 
‘hot -hown) along the west coast and over Oklahoma. 
Both sre shown in the mean circulation of figure 6A for 
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FicuRE 7.—Mean 700-mb. contours (solid) and height departures 
from normal (dotted), both in tens of feet for July 1-25, 1960. 
This represents the mean circulation prior to the formation of 
tropical storm Brenda (dotted track, east coast) which traveled 
up the mean trough (heavy vertical line). Also depicted are 
paths of hurricanes Abby and Celeste. Large dots indicate 1200 
GMT positions of storms, and numbers the July date of each 
position. 
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Figure 8.—Total precipitation (inches) for July 1960. Record 
monthly amounts for July were observed along the east coast. 
(From [2].) 


July 21-25. The Oklahoma trough was the stronger of 
the two, and heavy rain occurred in its vicinity during the 
23d and 24th. Subsequently a portion of this trough 
advanced to join the mean trough in which Brenda was 
spawned, while the trough near the west coast strengthened 
and advanced enough to bring showers and local tem- 
porary relief from critical burning conditions in forested 
areas of the Far West. 

Practically no rain had fallen over the western tier of 
States until this time. 
adjoining States and from Montana to Arizona (see fig. 8), 


Dry weather also prevailed over 


where the general accumulation was less than half an inch. 
Sheridan, Wyo., had 3.51 inches during the period from 
January to July, the driest like period ever recorded; and 
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Phoenix, Ariz., reported no measurable rain for 143 con- 
secutive days ending July 23. The dryness was a natural 
consequence of the persistent anticyclonic circulation over 
western United States. Migratory cyclones were com- 
pletely absent from the area west of the Continental 
Divide (fig. 5B). In addition, orographic precipitation 
was far below normal in the Pacific Northwest where the 


mean “jet”? was diverted far to the north (fig. 3). 


5. OTHER STORMS IN THE TROPICS 


Hurricane Abby was discovered near the Windward 
Islands. It reached hurricane intensity on the 10th, 
traveled steadily westward near 15° N. latitude, and 
skirted the northern coast of Honduras during the 14th. 
After moving inland on the 15th, the storm apparently 
dissipated and dropped from sight. On the basis of con- 
tinuity, however, it is likely that remnants of the storm 
were instrumental in the formation of hurricane Celeste 
on the 19th (see dashed continuation of Abby, fig. 7). 
Celeste moved west-northwestward well off the coast of 
Mexico, turned westward south of Baja California, on the 
22d, weakened to a tropical storm the following day, and 
dissipated thereafter. 

In the western Pacific two storms attained typhoon 


intensity. Typhoon Polly formed east of Luzon, P.I., on 
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July 19 and took a slow northerly track through the East 
China Sea and the Yellow Sea, reaching the coast of 
Manchuria on the 29th. Shirley also formed east of the 
Philippines, but took a more westerly course and struck 
Formosa on the last day of the month. 


REFERENCES 


1. L. M. LaMois, ‘‘Weather and Forest Fires,’’ Weekly Weather and 
Crop Bulletin, National Summary, vol. XLVII, No. 33, August 
15, 1960, p. 7. 

2. U.S. Weather Bureau, Weekly Weather and Crop Bulletin, Na- 
tional Summary, vol. XLVII, Nos. 27-31, July 4, 11, 18, 24, 
August 1, 1960. 

3. C. M. Woffinden, “The Weather and Circulation of June, 1960— 
A Hot Dry Month in the Southwest,’’ Monthly Weather Review, 
vol. 88, No. 6, June 1960, pp. 229-234. 

4. J. Namias, ‘“‘The Annual Course of Month-to-Month Persistence 
in Climatic Anomalies,” Bulletin of the American Meteorologi- 
cal Society, vol. 33, No. 7, Sept. 1952, pp. 279-285 (and an 
unpublished extension through 1954). 

5. W. H. Klein, “The Weather and Circulation of June 1957— 
Including an Analysis of Hurricane Audrey in Relation to the 
Mean Circulation,” Monthly Weather Review, vol. 85, No. 6, 
June 1957, pp. 208-220. 

6. C. 5S. Ramage, ‘‘Long-Period Circulation Anomalies and Tropical 
Storms,” Journal of Meteorology, vol. 17, No. 3, June 1960, 
p. 375, and Reply by E. M. Ballenzweig, pp. 375-376. 


Publications by Weather Bureau Authors 


J. K. Angell, “An Analysis of Operational 300 Mb. Transosonde 
Flights from Japan in 1957-58,” Journal of Meteorology, vol. 17, 
No. 1, Feb. 1960, pp. 20-35. 

K. Butson, ‘Florida Winter Weather: 1957-58,’ Weatherwise, vol. 
11, No. 2, Apr. 1958, pp. 58-59. 

N. L. Canfield, ‘‘Review of Atlas of 300-mb. Wind Characteristics for 
the Northern Hemisphere, by Lahey et al., Univ. of Wisconsin 
Press, Madison, 1960,” Bulletin of the American Meteorological 
Society, vol. 41, No. 7, July 1960, p. 399. 

C. B. Carney, “Freeze Damage to Plants,’ Weekly Weather and 
Crop Bulletin, National Summary, vol. XLVI, No. 21, May 25, 
1959, pp. 7-8. 

G. P. Cressman, ‘Hemispheric Nondivergent Barotropic Fore- 
easting,’’ The Atmosphere and the Sea in Motion, Scientifie Con- 
tributions to the Rossby Memorial Volume, Bert Bolin, Ed., 
The Rockefeller Institute Press in with Oxford 
University Press, New York, 1959, pp. 475-485. 

R. A. Dightman, ‘‘An Application of Climatology to Alfalfa Har- 
vesting,’ Weekly Weather and Crop Bulletin, National Summary, 
vol. XLVI, No. 28, July 13, 1959, p. 8. 

S. Fritz, “Satellite Meteorology,” in ‘United 
Report 1957-1960 Twelfth General Assembly, 
Union of Geodesy and Geophysies,’”’ Transactions of the American 
Geophysical Union, vol. 41, No. 2, June 1960, pp. 217-220. 

J. G. Georg, “Cranberry Crop and Weather in Wisconsin, 1959,” 
Weekly Weather and Crop Bulletin, National Summary, vol. 
XLVI, No. 41, Oct. 12, 1959, pp. 7-8. 

D. L. Harris, ‘Storm Surges,’ in ‘‘United States National Report 
1957-1960, Twelfth General Assembly International Union of 
Geodesy and Geophysics,’ Transactions of the American Ceo- 
physical Union, vol. 41, No. 2, June 1960, pp. 266-269. 

H. E. Landsberg, ‘Do Tropical Storms Play a Role in the Water 
Balance of the Northern Hemisphere?” Journal of Geophysical 
Research, vol. 41, No. 4, Apr. 1960, pp. 1305-1307. 


association 


States National 


International 


H. E. Landsberg, ‘Review of General Climatology, by H. J. Criteh- 
field, Prentice-Hall, Englewood Cliffs, N.J., 1960,’’ Science, vol. 
131, No. 3417, 24 June 1960, p. 1882. 

A. R. Long, “Dry Weather in Alabama, October 1958—August 1959,” 
Weekly Weather and Crop Bulletin, National Summary, vol. 
XLVI, No. 35, Aug. 31, 1959, p. 8. 

L. Machta, ‘“‘Meteorological Factors and Fallout Distribution,” 
pp. 33-56 in Low-Level Irradiation, Austin M. Brues, Editor, 
Publication No. 59 of the American Association for the Advance- 
ment of Science, Washington, D.C., 1959, 148 pp. 

L. Machta and R. J. List, ‘‘The Global Pattern of Fallout,’’ Chapter 
2 (pp. 26-36) in Fallout, A Study of Superbombs, Strontium 90 and 
Survival, John M. Fowler, ed., Basie Books, Ine., New York, 
1960, 235 pp. 

F. Ostapoff, “‘“Antaretic Oceanography” in ‘United States National 
Report, 1957-1960, Twelfth Assembly International 
Union of Geodesy and Geophysies,’”’ Transactions American ( e0- 
physical Union, vol. 41, No. 2, June 1960, pp. 257-258. 

J. C. Purvis, ‘“‘Weather and Peaches in South Carolina,’ Weekly 
Weather and Crop Bulletin, National Summary, vol. XLVI, No. 
2, Jan. 12, 1959, pp. 6-8. 

A. D. Robb, ‘“‘Rose Culture and the Weather,’ Weekly, Weather and 
Crop Bulletin, National Summary, vol. XLVII, No. 2, Jan. 11, 
1960, p. 8. 

L. H. Seamon, “‘Weather of the Year 1958,’’ Weekly Weathe, and 
Crop Bulletin, National Summary, vol. XLVI, No. 4, Jan. 26, 
1959, pp. 7-8. 

S. Teweles and F. G. Finger, ‘‘“Reducticn of Diurnal Variation ix the 
Reported Temperatures and Heights of Stratospheric Cons‘ ant- 
Pressure Surfaces,’’ Journal of Meteorology, vol. 17, No. 2, Apr. 
1960, pp. 177-194. 

H. Wexler, “Satellites et Meteorologie,”’ Met-Mar, Bu'‘etin 
Trimestriel de Liason et d’Information, No. 28, Direction «‘¢ la 
Meteorologie Nationale, Paris, April 1960, pp. 13-19. 


General 


U.S. GOVERNMENT PRINTING OFFICE: 1960 











